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Abstract. Various exclusive and inclusive decays of B mesons are studied, at present, with dedi- 
\ cated experiments at "B factories". In order to compete with the experimental accuracy, we need 

a reliable theoretical framework to compute strong interaction effects in a hadronic environment. 
I discuss how the separation of (perturbatively calculable) short-distance QCD effects from (non- 



perturbative) long-distance phenomena helps to obtain precise theoretical predictions. 

lnucl-th/03100081 CERN-TH/2003-236] 



=5 



INTRODUCTION 



> 
00 

O 
O 

The standard model (SM) for electroweak and strong interactions of elementary particles 
£T) , allows for precise theoretical predictions that are in remarkable agreement with present 

experimental results from high-energy experiments. However, an explanation of the 
SM and its parameters within a more fundamental theoretical framework remains an 
unresolved puzzle. Many questions in this respect (e.g. a grand unification of forces and 
matter, the origin of neutrino masses and its relation to the quark sector, the origin of CP 
violation, etc.) are related to the flavour sector of the SM, i.e. to the masses and coupling 
constants of the different quark and lepton species. 5-meson decays are particularly 
useful to explore the parameters of the flavour sector in the SM or possible extensions 
^ of it. The b quark has a relatively long lifetime, and provides many rare decay modes 

that have small branching ratios and are sensitive to details of the interactions at small 
distances. In this way one hopes to reveal indirect effects from physics beyond the SM 
even before the direct detection of new particles. 

From the QCD point of view, b quarks are interesting because they are the heaviest 
quarks that build pronounced hadronic bound states. The fact that the Z?-quark mass 
(mi, ~ 5 GeV) is large with respect to typical hadronic scales leads to new approximate 
symmetries that can be observed in the ^-hadron spectrum and decays. Furthermore, 
the heavy-quark mass provides a scale at which the strong coupling constant is still 
small, and perturbative computations of short-distance effects are possible. The short- 
distance dynamics has to be separated from (non-perturbative) long-distance physics 
related to the intrinsic QCD scale Aqcd- Technically this is achieved by an operator 
product expansion and from an effective field theory approach, respectively. In this way, 
5-meson decays can also be used to test and improve our theoretical understanding of 
QCD in different dynamical regimes. 
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FIGURE 1. Weak flavour transitions in the standard model. 



QCD AT THE ELECTROWEAK SCALE 

In the SM, transitions between different quark flavours are mediated by charged W 
bosons. The relative strengths Vij for flavour transitions g ; — > qj define the unitary 
Cabbibo-Kobayashi-Maskawa (CKM) matrix, which can be parametrized in terms of 
3 real angles and one CP-violating phase. On the other hand, flavour-changing neutral 
currents (FCNCs) can only be induced via loop diagrams (box or penguin topologies, 
see Fig. [l]). As a consequence, FCNCs are sensitive to the properties of virtual heavy 
particles in the loops, e.g. the top quark in the SM, or new particles in extensions of it. 



Effective Hamiltonian 

As in the case of the muon decay, we may "integrate out" the heavy particles 
(W, Z bosons, top quark, new physics) to arrive at an effective Hamiltonian (for a specific 
flavour transition), which has the schematic form 

#eff °= Gf 2^Q(aO ' ^ + terms suppressed by (1) 
t 

Here Gf is the Fermi constant, and C ; (/l) are effective coupling constants (Wilson 
coefficients) that encode the short-distance dynamics from physics above the scale 
/l = 0(m\y). The dynamics of the remaining five quark flavours, leptons, gluons and 
photons at energy scales below /i is described by a set of operators G\ (four-fermion 
operators, and operators that couple fermions to the chromomagnetic or electromagnetic 
field strength). For more details, see the review in 01. 



Matching 

The Wilson coefficients can be calculated by "matching" scattering amplitudes in the 
SM and the effective theory. This requires computation of QCD (and QED) radiative 
corrections in both cases; see Fig. El Since the strong coupling constant is small at the 
matching scale, the Wilson coefficients have a perturbative expansion, 

(o) a s {m w ) (i) 
Ci(m w ) =c)'+ c) ' + ... (2) 




effective theory: 




FIGURE 2. Example for radiative correction in the full and in the effective theory. 



At present, the SM matching coefficients are known at order aj (i.e. at two loops). 
Different models for physics at and above the electroweak scale yield different matching 
coefficients: 



» _ » 



c\' l> (m t ,m w , ...) +c ; -^(new physics 



(3) 



Therefore, experimental measurements of Wilson coefficients in weak decays test the 
SM and/or constrain new physics models. 



Resummation of large logarithms 

Radiative corrections to the matrix elements of the effective operators G\ in general 
would involve large logarithms lnm^//!, when calculated at a scale \x ~ my/- Since 
higher orders in perturbation theory also lead to higher powers of logarithms, the conver- 
gence of the perturbative series would be poor. This can be avoided by exploiting the fact 
that Wilson coefficients in the effective theory obey a renormalization-group equation, 
which allows the evolution of Q(m^) to C\{m^). At the low scale m^, the large loga- 
rithms (lnm^/ mw) n are explicitly resummed in C, (m^), and matrix elements of G{ only 
contain dynamics from energy scales below \x = nit,. In order to derive the evolution 
equations, one has to calculate the anomalous-dimension matrix jij mat describes the 
scale dependence and mixing of operators in the effective theory. Currently, the calcula- 
tion of three-loop anomalous dimensions is being completed (for a recent contribution, 
see Jl). 

Example: B — > X s y. A prominent phenomenological example, where the theoreti- 
cal machinery of the effective-Hamiltonian approach is relevant, is the inclusive rare 
radiative decay B — ► X s y. It provides a stringent test of the contributions to the Wilson 
coefficient Cj, which is related to the effective b — > sy vertex. The comparison between 
experiment and theory reads 

Exp.: BR[B^X s y] = (3.34±0.38) x 1(T 4 Jl (4) 
SM: BR[B^X s y] = (3.70±0.30) x 1(T 4 HH (5) 

The theoretical uncertainty is dominated by the renormalization-scheme dependence 
induced by the charm quark mass. The experimental and theoretical accuracy is already 
sufficient to put strong constraints on many new physics models (for a recent review on 
inclusive rare B decays, see |0]). 



QCD AND HEAVY-QUARK EXPANSION 

Because quarks are confined, experiments can only probe weak interactions in a 
hadronic environment. Technically, this amounts to considering hadronic matrix el- 
ements of the effective Hamiltonian <[T]). By construction, these matrix elements are 
sensitive to long-distance QCD dynamics, which is not accessible in perturbation theory. 
Nevertheless, some simplifications arise from the fact that the b quark mass is large 
compared to Aqcd- On the one hand, the strong coupling constant is small, OC s ( m b) <S 1, 
which implies that the dynamics at distances of order l/m/, is still perturbative. On the 
other hand Aqcd /nib "C 1 provides a small expansion parameter, and in the heavy-quark 
limit (nib °°) the number of independent unknown hadronic quantities may be less 
than in the general case. 



Heavy quark effective theory (HQET) 

The above observations can be formalized in terms of an effective theory (HQET) 
where - to first approximation - heavy b and c quarks are replaced by static 
colour/flavour sources, moving with a fixed velocity v^, 

h v (x)=e" n ^ l -^Q(x). (6) 

The first few terms in the effective Lagrangian, resulting from integrating out "small" 
spinor components and "hard" quark and gluon modes (with virtualities of order m^), 
reads 

^HQET = ^|/v-D+^+C m ( i u)^CT MV G^ V + ...|/ I v, (7) 

(where v • D — 0). The first term in this expansion is independent of the heavy-quark 
mass and diagonal in the heavy-quark spin. As a consequence, two new symmetries 
arise in the heavy-quark limit: heavy-flavour symmetry reflects the fact that the soft 
interactions of the b or c quark become the same, once the center-of-mass motion of 
the heavy quark is subtracted. Heavy-quark spin symmetry is related to the fact that soft 
interactions do not change the spin of the heavy quark. This has important implications 
for phenomenological observables such as the heavy hadron mass spectrum, or heavy 
meson transition form factors, to be discussed below. 1 The symmetries are broken by 
the subleading (kinetic and chromomagnetic) terms in the Lagrangian Q, as well as by 
perturbative matching coefficients for decay currents. For more details and references to 
the original literature, see for instance the review |0]. 



HQET also applies to heavy-to-light decays as long as the energy transfer to the light quarks and gluons 
is small. In this kinematic regime, exclusive heavy-to-light processes are described by an effective low- 
energy Lagrangian for heavy and light mesons that combines HQET and chiral perturbation theory 01 . 



Example: Extracting \V c b\ and \V u b\. A well-known application of the heavy-quark 
mass expansion and HQET is the determination of the CKM elements \V c b\ and \V u b\ 
from b — > c and b — > u decays. At leading order in the 1 /m-b expansion, the inclusive 
rates are just given by the partonic subprocess, which is calculable in perturbation theory 
(this has also been exploited in obtaining the theoretical prediction for B — > X s y in 
©). Power corrections to the inclusive B — > X rates start only at order l/m|, and the 
corresponding non-perturbative parameters can be fitted to moments of experimentally 
measured inclusive decay spectra. As a result, one extracts [9] 

I V^linci. = 0.0421 ±0.0013 , 

\Vub\mci. = 0.00426 ±0.00013 ±0.00050. (8) 

The theoretical accuracy is limited by the control on higher-order corrections in pertur- 
bation theory and in the heavy-quark expansion. For the extraction of \V u b\ a subtle issue 
is how to separate the b — > utv spectrum from the b — > civ background. 



Soft-collinear effective theory (SCET) 
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FIGURE 3. Kinematics for the decay B — > Xy at large recoil energy. 

The heavy-quark expansion can also be systematically applied to cases where the 
heavy quark decays into energetic light quarks and/or gluons, see Fig. El Examples are 
exclusive B — > nn decays iflOll . or the endpoint spectrum in inclusive B — ► X s y decay 
P . In addition to the HQET field © and "soft" (i.e. low -energetic) light quark and 
gluon fields, the effective theory (SCET) for these cases includes "collinear" quark 
and gluon modes whose energy is proportional to the heavy-quark mass. In particular, 
collinear quarks in SCET are described in terms of "good" light-cone components t, , in 
terms of which the Lagrangian that describes collinear quarks interacting with collinear 
gluons reads 

^coii = l\in--D+ (ip ± - m) — (ip L + m )\*±S. (9) 

Here n^_ is a light-like vector, which is determined by the jet axis or the momentum of the 
outgoing hadron(s), and n^_ is another light-like vector with (n+ -n_) = 2. Contributions 
to heavy-to-light decays from the "bad" spinor components of collinear quarks, as 
well as from interactions between soft and collinear particles are suppressed by either 
1 1 nib or by CC S (see e.g. lfl2l 13]). In this way SCET provides an elegant alternative 



to understand factorization theorems for various QCD processes 111 411 . Although the 



physical motivation for separating short- and long-distance physics in HQET and SCET 
is similar, there are important differences to be noted: 

• SCET gives no constraints for the light hadron spectrum. 

• Interactions between energetic particles in the final state with soft spectators in the 
B meson are mediated by particles with virtualities of order /i 2 = /i^Aqcd* which 
corresponds to a second short-distance scale in addition to m\. 

• The effective Lagrangian and currents in SCET are non-local. 

Example: Forward-backward asymmetry zero in B — > K*£ + £~ . The suppression of 
the "bad" spinor components for light quarks in SCET and heavy quarks in HQET leads 
to the reduction of independent form factors for heavy-to-light transitions at large recoil 
|[l5ll . As a consequence, to first approximation, the position of the forward-backward 
asymmetry zero in the dec ay spect rum of the exclusive B — > K*£ + £~ can be predicted in 
a model-independent way Ilia. 0711 in terms of the SM Wilson coefficients Cj and Cg for 
b — » sy and b — ► s£ + £~ transitions. First-order radiative corrections have been calculated 
in lfl8ll and lead to the estimate q\ = 4.2 ± 0.6, where q\ refers to the invariant mass of 
the lepton pair at which the forward-backward asymmetry vanishes. 

QCD AND HADRONIC EFFECTS 

After separating and calculating the short-distance QCD effects, we are left with matrix 
elements that encode the long-distance dynamics of quarks and gluons within hadrons. 
Since QCD perturbation theory cannot be applied in this case, we need alternative 
strategies to determine these hadronic parameters and obtain theoretical predictions for 
physical observables: 

• extract hadronic parameters from one experimental observable, and insert them into 
the prediction for another; 

• estimate them from non-perturbative calculations within lattice QCD; 

• estimate them using QCD (light-cone) sum-rule techniques; 

• tune phenomenological models to a subset of hadronic observables, and estimate 
other observables within that model. 

In all these cases, the ultimate challenge is to give a reliable prediction of the system- 
atic uncertainties. For the first option this apparently is a straightforward task, but only 
practicable if there are enough independent processes where one and the same hadronic 
quantity enters (which is the case, for instance, for parton distribution functions). In 
the lattice approach one has to deal with several (simultaneous) extrapolations: from 
finite lattice spacing to the continuum limit (a — > 0), from the "quenched approxima- 
tion" {rif — 0) to dynamical quarks, from quarks fitting on a finite lattice to realistic 
light- and heavy-quark masses. Here, progress is to be expected from improved lattice 
actions, increasing computer power and optimizing numerical algorithms, as well as 
from understanding the systematic effects for dynamical quark simulations better (see 
also C. Davies' talk at this conference). QCD sum rules are based on parton-hadron 



TABLE 1. Estimates for B-meson decay constant and B®-B® mixing parameter from lattice QCD and 
QCD sum rules (numerical values taken from |2(J ). 
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Lattice QCD 
QCD sum rules 



203(27) MeV 
208(27) MeV 



238(31) MeV 
242(29) MeV 



1.34(12) 
1.67(23) 



1.00(3) 

~ 1 



1-18(4)(> 2 ) 



duality, which is applied to appropriate correlation functions, from which properties of 
the lowest contributing resonance can be extracted (for an introduction to QCD sum 
rules, see lfl9lo . The suppression of contributions from higher resonances is achieved by 
a Borel transformation, and by relating the hadronic and the partonic representation of 
the spectral function above some threshold sq. Among others, systematic uncertainties 
arise from varying the Borel parameter and the threshold parameter within a "stabil- 
ity window", and from power corrections parametrized in terms of quark and gluon 
condensates. Finally, concerning model estimates, the only way to quantify systematic 
uncertainties often is to compare the predictions of sufficiently many different models. 

In the following I give some examples of hadronic parameters relevant to B physics. 
More details on the numerical estimates for various quantities can be found in KOll . 

fs and Be- Important phenomenological quantities that characterize the B meson 
are the decay constant and the B q -B q mixing parameter. They are defined as 

(0\q W b\B(p)) = if BqPll ; {B° q \0(AB = 2)\B° q )) = *-B Bq {n)fl q m\ , (10) 

respectively, where €?(AB = 2) is a four-quark operator that induces bq «-> qb. Numerical 
estimates for these quantities are summarized in Table [T] (where we considered the 
renormalization-group invariant quantity Bg q for convenience). The ratio t, defined in 
that table plays an important role in the analysis of the CKM triangle. The asymmetric 
error in the lattice estimate arises from the different ways of performing the chiral 
extrapolation to realistic light-quark masses, with or without including "chiral logs". 
We also remark that, in principle, 5-meson and D-meson decay constants are related by 
HQET. However, at present, both the experimental accuracy in measuring and the 
theoretical understanding of l/m c corrections are insufficient. 

HQET parameters. In HQET the masses of the ground-state pseudoscalar (P) and 
vector (V) mesons to order 1/m.Q accuracy read 

{ my ) = me+AK) + 2^1 -A 1 K) + (1-5,K ) )A 2 K) j • 

In the heavy-quark limit, pseudoscalar and vector meson belong to the same spin- 
symmetry multiplet characterized by a residual mass term K{m q ) that only depends on 
the flavour of the light degrees of freedom. The parameters X\ and A2 are related to 



the kinetic energy and chromomagnetic operator in (Q, respectively. 2 The parameters in 
(fTTTl can be determined from combing information from spectroscopy, lattice, QCD sum 
rules, and experimental data on inclusive 5-meson decays (see Il20ll for a summary of 
quantitative results). Notice that the definition of the quantities on the right-hand side in 
(fTTTl is renormalization-scheme and -scale dependent. 3 



Isgur-Wise form factor. In the heavy-quark limit, as a consequence of heavy-quark 
flavour and spin symmetry, all B — > £>W form factors are given in terms of one Isgur- 
Wise function, which is normalized at zero recoil ll30ll . e.g. 



(D\h v ,rh v \B)oc£(vv')tr 



l+f'l+f 



i . 



(14) 



Perturbative QCD corrections to £(1) can be calculated in HQET. Non-perturbative 
corrections to (fl4l start only at order l/m^. The estimate thus obtained for the B — > D* 
form factor relevant to the analysis of \V c b\ is lEoll 



l + 5 1/m 2 + 5 1/m 3+... 



0.91 



H0.03 
-0.04 • 



The knowledge of the form factor enables one to extract the matrix element \V c b\ from 
exclusive decay modes. In practice, one also needs an estimate of the slope £'(1) to 
extrapolate the experimental data to zero recoil. In this way one obtains |V c £|exci. = 



We included a correction 5y arising from the mixing between different vector mesons obtained from 
light degrees of freedom with spin-parity / = l/2~ and f = 3/2~, respectively. The mixing angle 
between members of different spin-symmetry multiplets vanishes in the heavy-quark limit, but may be of 
the order of 10° for the two lightest axial-vector D mesons 12111 . 

3 The recent measurement of new excited D and D s resonances 1 22] has renewed the interest in theoretical 
predictions for heavy-meson spectra 1 23]. The analogous formula to f° r the lowest-lying heavy scalar 
(S), and axial-vector (A) states reads 

I m A J u \ m 2m Q \ -X[{m q ) + {\-8 A {m q )))J l {mq) J 

In general, the HQET parameters for the two different (would-be) spin multiplets in dl It and dl 2t i 
are independent of each other. Additional constraints are obtained if one Taylor-expands all quantities 
around the chirally symmetric limit (i.e. vanishing quark condensate, (qq) — » 0), and if one assumes that 
higher-order terms in that expansion are suppressed ll24l l25ll . In this case the mass splitting between 
chiral multiplets is expected to be A'(m q ) — A(m q ) <x \(qq)\. The empirical values for this quantity 
are in fair agreement with phenomenological models based on spontaneous chiral symmetry breaking 
j26ll27ll28ll29ll . Furthermore, neglecting contributions proportional to (qq) or m q in terms that are already 
suppressed by 1 / mg, one expects 

X\(m q ) ~ k[(m q ) ee X\ , Xz{m q ) ~ k^triq) = X 2 , 5v(m q ) ~ 8 A (m q ) = 8 . (13) 

In other words, to first approximation 1 / mg corrections should not depend on the flavour and parity of the 
light degrees of freedom, which appears to be in good agreement with experimental data. Nevertheless, 
one has to keep in mind that neither l/m c nor (qq) are very small on hadronic scales, and higher-order 
terms in dl 1H13> may be non-negligible. 



0.0402 ± 0.0020 |9], which is in good agreement with the number obtained from in- 
clusive modes ©. 

Heavy-to-light decays. In order to extract information on heavy-to-light decays from 
exclusive channels, one needs (among others) the corresponding transition form factors, 
for instance to obtain C7 from B — > K*y or \V u t,\ from B — > n(p)£v. In contrast to 
the heavy-to-heavy case, the values of heavy-to-light form factors themselves are not 
restricted by a simple symmetry principle. Often, one uses a simple parametrization, 
e.g. for the relevant form factor for B — > ntv decays: 

fB-+*(J) = /+^*(°) (l5) 

J+ Kq) {l-q 2 /m 2 B *)(l-aq 2 /ml) ' K } 

where the first term in the denominator reflects the fact that for q 2 ~ m\ the form factor 
is dominated by the nearest vector-meson pole, and f+~**£0) — 0.2-0.3 and a = 0.3-0.6 
are fitted to lattice, QCD sum rules, or model estimates Il20ll . With increasing experimen- 
tal data on the exclusive decay spectra, which would also help to test and refine different 
model estimates, significant improvement is expected for the determination of \V u h \ from 
B — > n£v and B — > p£v decays. Currently, the result for \V u b\ from exclusive decays is 
somewhat smaller than the result from the inclusive analysis (but compatible within the 
rather large uncertainties) [9]. 

Other exclusive B decays into light mesons require additional non-perturbative input, 
as soon as radiative QCD corrections are considered. The first systematic treatment of 
this issue has been discussed in II 1 011 for the case of non-leptonic B — > %% decays. It has 
been shown that (to leading order in the 1 /mj expansion and including first-order a s cor- 
rections) the decay amplitude can be factorized into perturbative short-distance coeffi- 
cient functions and non-perturbative matrix elements that define heavy-to-light transition 
form factors, or light-cone distribution amplitudes for B mesons and pions. Comparison 
with experimental data from B — > %% and B — > TtK decays gives additional constraints 
on the CKM triangle (with some uncertainties coming from the parametrization of the 
5-meson distribution amplitude and from estimating I /mi, corrections). 

SUMMARY 

The different energy scales relevant to 5-meson decays imply that different dynami- 
cal aspects of QCD are probed: at the electroweak scale amplitudes are sensitive to the 
parameters that describe flavour transitions in the Standard Model or its possible ex- 
tensions. Perturbative QCD corrections are included by matching the full theory onto 
an effective Hamiltonian and evolving the corresponding Wilson coefficients to scales 
of the order of the heavy-quark mass. At this scale, strong interactions are still per- 
turbative. Furthermore, hadronic matrix elements can be expanded in inverse powers 
of the heavy-quark mass. This is formally described in terms of heavy quark effective 
theory (for b decays into charm quarks or low-energetic light quarks and gluons), of soft- 
collinear effective theory (for decays into energetic light quarks and gluons). Finally, the 
non-perturbative dynamics related to the intrinsic QCD scale is described in terms of 



hadronic parameters that have to be extracted from experiment or estimated from lattice, 
QCD sum rules or models. Reducing the theoretical uncertainties from each of the above 
dynamical regimes, and comparing the results with experimental data from B factories, 
enables us to test the flavour sector of the Standard Model, to find or constrain indi- 
rect new physics contributions, and to improve our understanding of perturbative and 
non-perturbative QCD. 



ACKNOWLEDGMENTS 



I would like to thank M. Nowak and A. Polosa for interesting discussions about the D s 
resonances, and T. Hurth for a careful reading of the manuscript and helpful comments. 



REFERENCES 

1. Buchalla, G., Buras, A. J., and Lautenbacher, M. E., Rev. Mod. Phys., 68, 1 125-1 144 (1996). 

2. Gambino, P., Gorbahn, M., and Haisch, U., hep-ph/0306079 (2003). 

3. Jessop,C, SLAC-PUB-9610(2002). 

4. Gambino, P., and Misiak, M., Nucl. Phys., B611, 338-366 (2001). 

5. Buras, A. J., Czarnecki, A., Misiak, M., and Urban, J., Nucl. Phys., B631, 219-238 (2002). 

6. Hurth, T., hep-ph/02 12304 (2003). 

7. Wise, M. B., Phys. Rev., D45, 2188-2191 (1992); Burdman, G, and Donoghue, J. R, Phys. Lett., 
B280, 287-291 (1992); Yan, T.-M., et al., Phys. Rev., D46, 1 148-1 164 (1992). 

8. Neubert, M., Phys. Rep., 245, 259-396 (1994). 

9. Schubert, K. R., talk presented at Lepton Photon 2003, Fermilab, 11-16 August 2003 

Jhttp : / /conferences . fnal . qov/lp2003/proqram/S6/schubert_s0 6 .pdf| ] . 

10. Beneke, M., Buchalla, G, Neubert, M., and Sachrajda, C. T., Phys. Rev. Lett., 83, 1914-1917 (1999). 

11. Bauer, C. W., Fleming, S., and Luke, M. E., Phys. Rev., D63, 014006 (2001). 

12. Bauer, C. W., Fleming, S., Pirjol, D., and Stewart, I. W., Phys. Rev., D63, 1 14020 (2001). 

13. Beneke, M., Chapovsky, A. P., Diehl, M., and Feldmann, T., Nucl. Phys., B643, 431-476 (2002). 

14. Bauer, C. W., et al., Phys. Rev., D66, 014017 (2002). 

15. Charles, J., Le Yaouanc, A., Oliver, L., Pene, O., and Raynal, J. C, Phys. Rev., D60, 014001 (1999). 

16. Burdman, G, Phys. Rev., D57, 4254-4257 (1998). 

17. Ali, A., Ball, P., Handoko, L. T., and Hiller, G, Phys. Rev., D61, 074024 (2000). 

18. Beneke, M., Feldmann, T., and Seide l, P., Nucl. Phys., B612, 25-58 (2001). 

19. Colangelo, P., and Khodjamirian, A., hep-ph/0010175 (2000). 

20. Battaglia, M., et al., hep-ph/0304132 (2003). 

21. Kilian, U., Korner, J. G, and Pirjol, D., Phys. Lett, B288, 360-366 (1992); Di Pierro, M., and 
Eichten, E., Phys. Rev, D64, 114004 (2001); Abe, K., et al., hep-ex/0307021 (2003); Cheng, H.-Y, 
hep-ph/0307168(2003). 

22. Aubert, B., et al., Phys. Rev. Lett, 90, 242001 (2003); Besson, D., et al., hep-ex/0305017 (2003); 
Abe, K., et al., hep-ex/0307041 (2003). 

23. See also the numerous contributions to the proceedings of this conference. 

24. Bardeen, W. A., Eichten, E. J., and Hill, C. T., hep-ph/0305049 (2003). 

25. Nowak, M. A., Rho, M., and Zahed, I., hep-ph/0307102 (2003). 

26. Nowak, M. A., Rho, M., and Zahed, I., Phys. Rev, D48, 4370-4374 (1993). 

27. Bardeen, W. A., and Hill, C. T., Phys. Rev, D49, 409-425 (1994). 

28. Ebert, D., et al., Nucl. Phys., B434, 619-646 (1995), Phys. Lett, B388, 154-160 (1996). 

29. Deandrea, A., et al., Phys. Rev., D58, 034004 (1998), and hep-ph/0307069 

30. Isgur, N., and Wise, M. B., Phys. Lett, B232, 113 (1989) and B237, 527 (1990). 



